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Abstract 

Recently we reported some interesting features of the Wolff's algorithm 
behavior when applied to the site-bond-correlated Ising model. Our main 
results were that a stronger correlation diminishes the autocorrelation 
time but it does not change or affect the critical dynamic exponent. In 
this paper, we analyse the Wolff's cluster structure and how it varies with 
the spatial correlation. The fractal dimensions are determined for several 
values of the magnetic atoms concentration and spacial correlation, giving 
an estimative of the critical dynamic exponent. 
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With the advent of the cluster algorithms, a great improvement in the simu- 
lations of magnetic spin systems has been possible (Swendsen 1987, Wolff 1989). 
When applied to the usual dilute Ising model (Hennecke and Heyken 1993), such 
algorithms exhibit a performance which is better than those of single-spin flip 
techniques. It was verified that as the concentration of non-magnetic atoms in- 
creases the corresponding autocorrelation time decreases. Recently, we obtained 
very similar results about the performance of Wolff algorithm in the site-bond 
correlated Ising model (Campos and Onody 1997, Campos and Onody 1998). 
It was shown that the algorithm becomes even more robust in the presence of 
spatial correlation. This a very important feature, since for the single-spin flip 
algorithms (Metropolis, Glauber) the correlation weakens their performance, 
and so a poor statistical data can result. In this letter we deal with the analysis 
of the Wolff's clusters structures and their importance in obtaining the critical 
dynamic exponent. 

The site-bond correlated Ising model (hereafter, SBC model) was proposed in 
a way to explain additional effects not predicted by the simple dilute Ising model, 
and experimentally verified in the antiferromagnetic compound KNi^Mgi^xPs 
(Aguiar, Engelsberg and Guggenheim 1986). The main attribute of this model 
is the existence of a local spatial correlation. In the SBC model, the presence of 
impurities in the neighbourhood of a given (nearest neighbors) pair of magnetic 
atoms can modify the strength of the exchange coupling constant between these 
two atoms. Moreover, in the limit of maximum correlation a, the exchange 
interaction can be even suppressed. 

The SBC model Hamiltonian is defined as follows (Aguiar, Moreira and 
Engelsberg 1986) 



i,S 

where aj = ±1 and 6 denotes an elementary lattice vector. The exchange 
interaction Ji^i+s is given by 
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J i,i-\-5 — J^'. 
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where J > 0. The random variables Si can take the following values: one with 
probability C and zero with probability 1 — C, where C is the concentration of 
magnetic atoms. The parameter a correlates the interaction between sites i and 
i + 5 with the magnetic occupancy of the sites i — d and i + 26. The uncorrelated 
dilute Ising model is re-obtained in the limit a = 1. For < a < 1, the bond 
between i and i + 5 is only weakened by the absence of a magnetic atom at 
i — 5 OT i + 2S. The limit a = corresponds to the maximum correlation, 
i.e., two magnetic first neighbor sites are connected by an active bond only if 
their nearest-neighbor sites along the lino joining them are also present. This 
maximum limit (a = 0) defines at T==0 a new percolation problem, with the 
concentration threshold pc = 0.7405 (Campos, Pessoa and Moreira 1997). For 
a 7^ the critical concentration is like that the usual site percolation. 

Following the Fortuin-Kasteleyn (Fortuin and Kasteleyn 1969) arguments 
and the definition of the model hamiltonian, we can derive the bond activation 
probabilities required in the construction of the clusters. The probabilities are 
defined as follows (Campos and Onody 1997): 

Pi.i+5 = li ai ^ ai+s 

Pi.i+s = 1 - e~'^^ if ai = ai+s and ei-sSi+25 = 1 (3) 
Pi,i+S = 1 - e~'^°'^ if ai = ai+s and Si-sSi+2S = 

between sites i and i + S. 

We simulated the SBC model for some values of concentration C and corre- 
lation a. Due to the large geometric fluctuations, an accuracy statistical study 
and a great computational efforts were performed. In Table 1, we show our 
results for the critical dynamical exponents z for some values of the parame- 
ters. These results were obtained by a log-log plot of the autocorrelation time 
T versus L, since as expected r obeys a power law like t L^. 

The values of the critical dynamical exponent z for some values of correlation a 
and concentration C. 





C = 0.80 


C = 0.70 


a = 0.1 


0.26 


0.10 


a = 0.3 


0.34 


0.11 


a = 0.5 


0.38 


0.33 


a= 1.0 


0.40 


0.38 
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As we can see, there is an increasing of the z exponent with augmented 
values of correlation. This behavior is verified for all concentration values. In 
order to compare the performance of the Wolff dynamics with others, we need 
to rescale r (Wolff 1989), since one Wolff step has a computational cost of the 
cluster size \S\. So, we define the Wolff autocorrelation time by the relation 

r -rM (5) 

where d is the lattice dimensionality, and CL'* is the total magnetic mass of the 
system. Here, we restrict our study to the case d= 2. Like r, the Wolff time 
also exhibits a power law of the form ^ L^"" , where Zyj is the Wolff critical 
dynamical exponent. 

Another way to estimate the vahics, consists in regarding a power law be- 
havior for the mean cluster size (151) ~ L^-^ , where Df is the fractal dimension 
of the Wolff clusters. In this case, we can calculate z^ through the equation 

z^=z-{d-Df). (6) 

In the case of the usual Ising model, the advantage of this procedure is the 
removal of one source of error, since {\S\) = (x) ~ , where x is the magnetic 
susceptibility and ^ is the ratio of the exactly known Ising exponents (Wolff 
1989). 

In a recent work, wc verified no dependence of the z^, vahics with the pa- 
rameters C and a (Campos and Onody 1998). We found a critical dynamical 
exponent value around zero. This result indicates that the Wolff algorithm 
circumvents almost completely the critical slowing down phenomenon, and so 
making possible a more accuracy statistical analysis, since there is no statistical 
correlation between the sucessive configurations generated by the algorithm. 

But this an intriguing thing. With this assumption, the scenario presented 
for the z exponent is not valid for z^. The only way to support these results 
regards the changes occurred in the fractal dimension of the Wolff clusters with 
the parameter a in order to keep constants the Zw values . The study of the 
cluster structures are our main focus on this work. 
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In Figure 1, we present two pictures of the Wolff clusters produced in different 
generations. On the top of the figure, we have a = 0.1 and C = 0.7. On the 
bottom, a = 1.0 and the same value for C. The configurations were generated 
in the respective critical temperature. As we can see, the clusters for a = 0.1 are 
more compact than those for a = 1.0. This means that the value of {d — Df) is 
smaller for a = 0.1 than for a = 1.0. For this reason, the tendency of increasing 
of z with a, is not predicted for the values. 
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